WRIGHT FUNCTIONS GOVERNED BY FRACTIONAL 
DIRECTIONAL DERIVATIVES AND FRACTIONAL ADVECTION 

DIFFUSION EQUATIONS 



MIRKO D'OVIDIO 



Abstract. Wc consider fractional directional derivatives and establish some 
connection with stable densities. Solutions to advcction equations involving 
fractional directional derivatives arc presented and some properties investi- 
gated. In particular wc obtain solutions written in terms of Wright functions 
by exploiting operational rules involving the shift operator. We also consider 
fractional advection diffusion equations involving fractional powers of the nega- 
tive Lapla<;e operator and directional derivatives of fractional order and discuss 
the probabilistic interpretations of solutions. 



1. Introduction 

In this work we study the solutions to fractional Cauchy problems involving the 
operator 

n 

(1-1) (a-V) = ^afea,, 

fc=i 

where V = {dx^, ■ ■ ■ ,dx^) and ||a|p = + • • • + a^^ = 1. In particular, we are 
interested in studying the solutions to the fractional advection equation 

^ + (a-V)"V = 0, a,/3e(0,l] 

for the scalar field tjj = ip{x,t) with constant velocity a. By exploiting operational 
rules involving the shift operator, we obtain some representations of the scalar field 
^ in terms of Wright functions. Prom this, we arrive at the fractional advection 
(convection) diffusion equation 

{jP in 

— ^-i-A/w-ia-S/rw, ^,a,p€{0,l] 
(where A = V • V is the Laplace operator) and the transport equation 

— w = -(-a • V)"w - A(/ - K)w 

where K is the Frobenius-Perron operator corresponding to some transformation. 
For the above equations we study the solutions from a probabilistic point of view. 

Before starting with an overview of the work we introduce some notations which 
turn out to be useful further in the text: 
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• u ~ u(x, t) is a general solution to a general boundary value problem, 

• M is the Laplace transform of w, 

• u is the Fourier transform of u, 

• a • X = X]fc=i '^kXk and ||x|p = x • x, 

• is the Dzerbayshan-Caputo time-fractional derivative, 

• 9f and 9" are the Riemann-Liouvillc time- and space-fractional derivatives, 

• dx = d/dx and dt = d/dt, 

• ha is the density law of the stable subordinator S)'^ , t > 0, 

• is the law of £f , i > which is the inverse to 9)^ . 

After some preliminaries and auxiliary results, we introduce, in Section 3, the 
rule 

(1.2) e^d^f{x)^ f{x + 

where the shift operator exp C,dx plays an important role in studying fractional 
powers of and (1.1)- By exploiting such a rule we define the fractional directional 
derivative 



(a • V)" = (l - e-^^-^-^)) d7(s) 



where 

d-f{s)/ds = — ^ 

1 (1 — a) 

is the Levy measure of a stable subordinator. In Section 4, we study fractional 
equations of the form 

(1.3) (^^-f (a.V)"^u„,^(x,t) = 0, (x,<)eM!;x(0,+^) 

for a, /3 e (0, 1]. In particular we show that 

Uq,^(x, i) =W|(a-x,i) 
is the solution to the equation (1.3) subject to the initial and boundary conditions 

Ua,p{x,0) = S{xi) X • • • X 5{Xn), WQ,^(0,t) = 

where 

(1.4) U?{x,t) ^E^'SiSj"^), xeR,t>0 

is the density law of the composition involving a stable subordinator j^" and an 
inverse process £f . If a = /3 g (0, 1), then the function (1.4) becomes the Lamperti's 
law (see for example James [16]; D'Ovidio [9]) and the solution to (1.3), for a — /3, 
can be explicitly written as 

sin/?7r (a-x)''-ii^ 
u/3,/3(x,i) - ^^(a.x)2/^ + 2(a.x)^i^cos/37r + t2/3' ^ ^ 

This is the case in which the density law of the composition of processes appearing 
in (1.4) coincides with the law of the ratio of two independent stable subordinators 
jSjf , j — 1, 2, that is 



(1.5) X iJof/2.Qf, t>0. 
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A special case of (1.3), for (3 S (0, 1), is the time-fractional equation 
(af + (a- V))t;;3(x,i) = 0, (x,i) eR^x (0,+oo) 
subject to the initial and boundary conditions 

t.^(x,0) -,5(x), vpiO,t)=t-'"' 
for which we show that the solution can be written as 



where 



2TTi 



Ha 



(1-6) -g fc!r(,fc + p) ' M>-i,peC,zeC 

(_ff a denotes the Hankel path) is the Wright function which has been introduced and 
investigated by the British mathematician E. Maithland Wright in a series of paper 
starting from 1933. By normalizing Vfj we obtain the n-dimensional probability law 
pp whose one-dimensional marginals coincide with the density law of the hitting 
time £f = inf{s > : i^f ^ (0,*)}, t > which is the inverse to the stable 
subordinator j^f , t > 0. In particular, we get that 

Pr{iX^ edxi,..., „xf e dxn} 
^^("'^^^ dx,---dx„ 

is the law of the process ,'^p{t) — {lX^ , . . . ,n^f ) which can be regarded as the 
inverse of 

J^^i(x) = ^ X li^f^ + . . . + ^ X „iof^^ 

in the sense that 

Pr{3^p{t) < x} = Pr{jr^(x) > t}. 

Due to the fact that i^f , i > 0, has non-negative increments, the multi-parameter 
process .^(x), x e M" , possesses non-decreasing paths. 
For the time-fractional equation 

(9f + (a • V))il^(x, i) = 0, /3 e (0, 1), (x, t) e M!^ x (0, +^) 

subject to the initial and boundary conditions 

ii;5(x,0) = (5(x), iX'^iO^t) neN 

we obtain that 

ii^(x,<) = r-"^-iw-_^,„_„^ (^lir) 

and, for n, m G N, we show that 

il^(x,i)*il^(y,t)=il^+"(x + y,t) 

where * stands for the Laplace convolution over t. Furthermore, for n = 1, we 
obtain the density law of £f , t > 0. 
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In Section 5, we present the solutions to the time- fractional problems involving 
the operator 

(71 n „2 
i=0 j=0 ^ J 

and, for a = 2, wc find out some connection with Gaussian laws. In particular, we 
show that the solution to 

^ - (a • V f\ .9(x, t) = 0, /3 e (0, 1), (x, t) e M" x (0, +c^) 

subject to the initial condition (7(x,0) — (5(x) is written in terms of the Wright 
function (1.6) as follows 

(1-9) ff(x,i) = ^M^_|.i_|(-^ 

Furthermore, formula (1.9) can be written as 

(1.10) g(x,t) = / -^^lp{s,t)ds 

where 

(1-11) hM^l-pW-f^.i-p{~^) 

is the law of the inverse process £f , t > 0. 

Finally, we arrive at the fractional advection diffusion equation. We show that, 
for ■& £ (0, 1), a,f3 £ (0, 1) and a e M" such that ||a|| = 1, the solution w = w{x,t) 
to the fractional equation 



for (x, t) £ M" X (0, +oo), is given by 

w(x,t)= / dz I (is72ij(x — s, z) ft,Q,(a • s, z) Z^(z,t) 

where 72i?(x, t) is the law of the isotropic stable Levy process S2ij(i), ha{x, t) is the 
law of the stable subordinator Sjf and lp{x,t) is the law of the inverse process £f . 
The distribution w represents the density law of a subordinated E"-valued stable 
process with stable subordinated drift given by 

W(t) -S2^(£f)+aio^,,, t>0. 
We also show that the stochastic sohition to the transport equation 

(1.12) ^=_(a.V)"w-A(/-if)w 

where a,/3 € (0, 1] and K — e^^'^ is the shift operator, is given by 

=N(£f)+ai32,. 

where Nt is a Poisson process. If a = 0, then K = B becomes the backward 
operator and therefore w becomes the probability of the fractional Poisson process 
N(£f ) (see for example Beghin and Orsingher [5] and the references therein). 
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2. Preliminaries and auxiliary results 
For TO — 1 < a < m, m e N, we define the Weyl's fractional derivative 



T{m — a) 



gm 

r(TO — a) 



{x - ' 



which is named Riemann-Liouville derivative if a; > whereas, for to — 1 < /3 < to, 
we define the Dzerbayshan-Caputo fractional derivative as 



it) 



^m—f3— 1 

T{m- (3) 



* dTf{t) 



1 



r(TO - Jo [t - 5)"+"^ 



t > 0. 



The symbol * stands for the Fourier convolution and for a = /3 = 1, both fractional 
derivatives become ordinary derivatives. Furthermore, we recall the following con- 
nection between the above fractional derivatives ([27]) 



(2.1) 



gl3 f ™-l 



k=0 



t=o+ r(fc- /3 + 1) 



which leads to the following relation between Laplace transforms 

, ^ m— 1 



dtp 



(A) = af /(A) - J2 9tfit) 



k=0 



t=0+ 



k-l 



where af /(A) = A^/(A). 

Our aim in this section is to introduce the operational solutions to the fractional 
equation 



(2.2) 



subject to the initial and boundary conditions Up{x,0) = (po{x) and Up{Q,t) ~ 
with a,/3 € (0, 1]. The problem is to find an explicit form for the convolution 



(2.3) 
where 

(2.4) 



li;^{x,t)^Ep{-tPd^)ipoix) 



27rz J Ha - z 



k=0 



r(^fc + i)' 



> 0, zeC 



{Ha is the Hankel path) is the Mittag-Leffler function. We first observe that u{z) = 
Ei3{~wP z)H(z), where u; > and H{z) is the Heaviside step function, is the 
fundamental solution to the fractional relaxation equation 



(2.5) 



{z)+w^u{z) = 0. 



For (/3o = 5, we get ^(C) = 1 and thus, being d^Uf^{^,t) = ^°'Up{£^,t) the Laplace 
transform of the Riemann-Liouville derivative of lA^ , we obtain the Laplace trans- 
form 



(2.6) 
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As we can immediately check, from (2.5) or the well-known fact that is an 
eigenfunction for the Dzerbayshan-Caputo derivative, we reobtain 

From (2.6) and the fact that (see for example Samko et al. [27]) 

EB(-z'^a°') = — / —5 —5 -, a > 0, 

' TT 7o + 2w/5 cos ^vr + 1 

we can write 

(2 7) UHi t) r e~^^^" dx 

If we formally rearrange (2.3) from (2.7) we get that 

which takes, for a = /?, the following form 

(2.9) u,{x,t) = —j^ '^%2,+2./^i/'cos/37r + t2^^"(^-^) 

by taking into account the rule (1.2). Thus, from (2.9), for ipi^ ~ 5, we arrive at 

P ' TT a;2/3 + 2xHP cos /Jtt + t^P 

which is the solution to the equation (2.2) or the inverse Laplace transform of 
(2.7) for a = (3. For a well-defined function (^g, the explicit representation of the 
convolution (2.3) is therefore given by (2.8) or, for a = /3, by (2.9). 

The solution Up{^,t), x > 0, t > for a, € (0, 1), can be regarded as the law 
of the composition io^^j, t > 0, where Sjf, t > is an a-stable subordinator (see 

Bertoin [6]) with Laplace transform 

(2.10) Eexp-^.^^ cxp-t^" 

and £f , t > 0, is the inverse to the /3-stable subordinator i^f (see for example 
Baeumer et al. [3]) for which 

(2.11) = Ep{-tPO- 
We also recall that 

(2.12) Eexpi^i^r = exp (-t(-i^)") = exp (-il^l^e"'"^ i 

is the characteristic function of a totally (positively) skewed stable process. The sta- 
ble subordinator Sjf is a Levy process with non-negative, independent and station- 
ary increments whereas, the inverse process £f has non-negative, non-independent 
and non-stationary increments as pointed out in Meerschaert and SchefHer [20] and 
turns out to be fundamental in studying the solutions to some time-fractional equa- 
tions. We refer to £" as the inverse of Sjf in the sense that P{2f < x} = P{^" > 
t}. We notice that, from the fact that Sjf has non negative increments, the process 
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£" — inf{s > : i^" ^ (0, t)} is an hitting time. Let us write ha for the law of i^" 
^ for the law of £f , 



and If^ for the law of £f . The governing equations are known to be 



(2.13) (^^^+d:^h^{x,t)^0, x>Q,t>0 

subject to the initial and boundary conditions 

ha{x, 0) = S{x) 
ha{0,t) = 

and 

(2.14) (^d^ + ^^lp{x,t)=0, x>0,t>0 
subject to 



(2.15) 



laix, 0) = S{x) 

/^(Oji) = r(i-;3) • 



In view of (2.1), formula (2.14) with conditions (2.15) can be rewritten in the 
compact form 

(2.16) (^_ + _j/^(2;,t)^0, x>0,t>0. 

From (2.5), by considering that dxlpi^jt) — £,li3{£,,t), we obtain that lj3{(,,t) — 
Ejsi—t^S,) which is the Laplace transform (2.11). We notice that (see D'Ovidio 
[10]) 

tha{t,x) xha(x,t) ^ (r\ ^\ + ^ n 

where la can be written in terms of the Wright function (1.6) as in (1.11). Thus, 
we can write the solution to (2.2) as 

(2.17) Up{x,t) = / ha{x, s)li3{s,t)ds ~ — / s /^(s, a:) Z^(s, t) ds 

Jo ^ Jo 

and, from (2.10) and (2.11), we can write 

/•oo 

Jo 

which is in accord with (2.6). 

From the discussion above we arrive at the following fact 



(2.18) Ep{~t^A")= dsUj^{s,t) 

Jo 



which holds for some well-defined operator A. In particular, for A — dx, from 

(2.18) we write (2.3). 

We also introduce the homogeneous distribution 

(2.19) zl = f H{z), z e M, > -1 

r(i + 77) 
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(where H is the Heaviside function) which is locally integrable in M \ {0}. The 
Fourier transform of (2.19) is written as 



/ + 00 
-OO 

and can be obtained by considering the fact that 

Furthermore, as a straightforward calculation shows 
(2-21) I e^'^zldz^^. 

3. Fractional Directional Derivatives 

Let {X, d, fi) be a metric space, that is {X, d) is a locally compact separable 
metric space and fi is a Radon measure supported on X. Let L^di) be the space 
of all real-valued functions on X with finite norm 

(3.1) \u{x)\Pdfi{x), l<p<cx). 

Let Tt be a strongly continuous semigroup on LP{^i) and A be the infinitesimal 
generator of Tt, that is 

Ttu - 



(3.2) lim 

t-i.O+ 



t 



Au 







for all u £ D{A) where 

Dom{A) — {u € Lpdj.) : the limit (3.2) exists}. 
We formally write 

(3.3) Tt = e*-^, t > 

and T" will denote the semigroup for the fractional generator A". 

The fractional power of operators has been examined in many papers since the 
beginning of the Sixties. A power a of a closed linear operator A can be represented 
by means of the Dunford integral 

(3.4) A"" =^ — [ dX X'' {X - Ar\ ma}>0 

2m Jr 

under the conditions 

(i) A G p{A) (the resolvent set of A) for all A > 0; 
(m) \\X[XI + A)-^\\ < M < oo for ah A > 

where T encircles the spectrum (j{A) counter-clockwise avoiding the negative real 
axis and A" takes the principal branch. The reader can find information on this 
topic in works by Balakrishnan [4], Hovel and Westphal [15], Komatsu [17], Kras- 
nosel'skii and Sobolevskii [18], Watanabe [30]. For U{a} G (0, 1), the integral (3.4) 
can be rewritten in the Bochner sense as follows 

(3.5) A" = ''^^^ r dXX^-\X + A)-^A. 
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From the proof of the HiUe-Yosida theorem we have that 

/"OO 



By inserting this expression into (3.5) we get that 



e 



where 

/•oo 

/ dss^'-^e-' ^r{a) 
Jo 

is the Gamma function for which the fohowing weU-known property holds true 

T{a)T{l~a) = 

sm na 

and 

poo 

(3.6) / ds s-"'e-'^ = r(l - a)A"-^ 
Jo 

which holds only if < a < 1 (the reader can also consults the interesting book by 
Renardy and Rogers [26]). The representation (3.5) can be therefore rewritten as 

A" - A'^-^A, a e (0, 1) 

by considering (3.6). After some algebra, by exploiting the same calculation, we 
can also arrive at the representation 

(3.7) A" ^A^-A"-", n-l<a<n, nGN. 
The semigroup Tf of the translation operator 

A^^ 
dx 

is strongly continuous in Lp{^) with diJ,{x)/dx = e~"^, w > and \\Tt\\p — e""*/^ 
(see for example [8; 27]). In what follows we will study the fractional power of the 
directional derivative operator A — (a • V). 

The directional derivative of fractional order a E (0,1) U (1,2] has been first 
introduced by Meerschaert et al. [21] in connection with models of anomalous diffu- 
sions. By using the Fourier transform convention the authors define the fractional 
derivative D^^ of a well defined function / : M" !—> M by requiring that 

(3.8) D^(0 = /(0 / {-t^-erM{de), ^ eM", ae (0,1)U(1,2] 

J||e||=i 

where M{d6) is a probability measure on the unit sphere 

§"-1 = {0 e M" : ||6»|| = 1}. 
In the one- dimensional case, 6 = ±1 so that \\d\\ 1. Formula (3.8) becomes 

with p + q = 1 and therefore the inverse Fourier transform yelds the governing 
operator of a one-dimensional asymmetric stable process (or symmetric if q = 1/2). 

In order to define fractional directional derivatives and study fractional equations 
involving the operator (1.1) we first introduce the rule 

(3.9) e<^^f{x)^f{x + 
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for an analytic function / : M i— > M and ^ G M. The operational rules (3.9) plays 
an important role in quantum mechanics ([2]) and formally, can be obtained by 
considering the Taylor expansion of the analytic function / near x written as follows 



A;=0 



Thus, we get that 



fi^ + C) = E ll ^^'^(^) = E It ^'/(^) = e<^-f{x). 
k=o ' k=a 

However (Feller [13]; Hille and Phillips [14]), the Taylor series and therefore the rule 

(3.9) , can be extended to the class of bounded continuous functions on (0, +oo). 
The reader can find further details on the Taylor expansion for generalized functions 
in works by Estrada and Kanwal [11]; Anastassiou [1] and the references therein. 

In our view, from (3.7), (3.6) and (3.3) the fractional power of ,4 can be written 
as follows 

(3.10) A" = AA"-'^ = A— r/ dss-^Ts 

r(l - a) Jo 

and therefore, based on the discussion made so far, we introduce the following 
representation of the fractional power of the directional derivative (1.1). 

Definition 1. For < a < 1, a e M" such that ||a|| = 1 and V = {dxi, ■ ■ ■ ,dx„) 
we define 

(3.11) (a. V)" := rdss"^e'^(-^\ 

r(l - a) Jo 

Remark 1. We notice that, after some algebra, from (3.11) we get 

(3.12) (a • V)" = {l- e-^^-*-^)) ^7(5) 
where 

d-f{s)/ds ^ — 

i (1 — a) 

is the Levy measure of a stable subordinator (Bertoin [6]). Indeed, a Berstein 
function b : (0, +00) i-> M and therefore such that & > and (— l)'^9^6(a;) < for 
all X > and fc G N, is a Berstein function if and only if 



(3.13) b{x) = ci+C2X+ (1 - e^^)dTO(s), x>0 

"'0 

for constants ci,C2 > and a non-negative measure m on [0,cxd) satisfying 

(3.14) / (s A l)dm(s) < 00. 



Furthermore, & is a Bernstein function if and only if there is a convolution semi- 
group, say (f)t, such that 



(3.15) / e-^"d0t(s) = e 



-tb{x) 

Here, we consider b{x) = x" with ci — C2 = and thus, m = 7 and 0t(s) = ha{s, t) 
is the law of io". 
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Remark 2. In the one-dimensional case, from (3.11), for / : (0, +c«) i— M, we 
obtain that the fractional derivative 



d_ 
dx 



1 



r(i - a) Jo 



where —dx is the generator of Tg, takes the form 



d^fix) 



d_ 
dx 
d_ 
dx 
d_ 
dx 



1 



1 

r(i-«) 
1 

r(i-«) 
1 

r(i - a) 
d 



ds s "/(x — s) 

c 

ds s~°'f{x — s) 
dsfis) 



a G (0,1) 



T{l-a)dxJo (x-s)"' 

which is the Riemann-Liouville fractional derivative on the half line (0, +oo) ([27]). 
For / : (— oo, +oo) !—> M, we also obtain that 

d 



d^fix) 



dx 
d_ 
dx 



1 



1 



1 r°° 

dss-"f{x-s) 

, ae(0,l) 



r(i - a 
d 



-V Vp(x,i) + cDV(x, t) 



ds (x — 

o 

dsfjs) 

r(l-a) dx (x-s)" 

which is the Weyl derivative on the whole real line (— cx), -(-cx)), see for example 
Samko et al. [27]. 

4. Fractional advection equations 

This section has been inspired by the works by Meerschaert et al. [21, 23] where 
the authors studied the advection-dispersion equation 

dp 
di 

c > 0, introduced to model anomalous dispersion in ground water flow. The frac- 
tional operator D^^ is defined in the space of Fourier transforms in the sense that 
D^^p(x,t) is given as inverse of (— • 9)"p{^,t) with \\9\\ — 1. Thus, the operator 
D^j generalizes the directional derivative D\,j = 6 W. Here we study the advection 
equations in which the fractional directional derivative has been previously defined 
in (3.11) and find results based on the operational rule (3.9). 

We will show that, for the operator A = (a • V), the solution to the space 
fractional equation 

(I 

subject to the initial condition 

Ua{x,0) = uo{x) 
can be represented as the convolution 
(4.2) Ua,{x,t)=Tt"uo{x) 



A' 
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where 

Jo 

The semigroup d(j)t{s) — dsh^is, t) represents the law of the positively skewed Levy 
process 55", t > Q and plays, in this context, the role of 7(5) in (3.12). Due to the 
fact that jo" t (which is the elementary subordinator) for a — >■ 1, we obtain 
that 

lini ha{x, s) = 5{x — s) 

a— ^1 

and therefore 

u{x,t) = TtUa{x) 
is the solution to (4.1) for a = 1. 

We start our analysis by noticing the following fact. 

Lemma 1. For the solution to 

(4.3) (^^^+d:^u^^Q, x>0,t>Q 

subject to the initial and boundary conditions 

Uq(.T,0) = uo(x), Ua{0,t)^0, 

we have that 

(4.4) Uaix, t) = e~*^° Uf){x) = Ee-^'*°^- uo{x). 
Proof. By considering the function u : (0, +00) 1— > M with 

(4.5) u\^^^ = 0, 
we get that 

/+00 
dxe^^-d^uix) 
-00 

r+oo 



dxe'^""—^— / dzz-°'e-''^-u(x). 

r(i-«)yo ^ ' 

From the Foubini's theorem and by taking into account (3.9), we obtain that 

1 rcc r+oo 

diuiO = Y(l-a) J„ / dxe'^^d,e-'^^u{x) 



1 



Til -a) 
1 



-00 
+00 



dzz " j dx e^^^ dx u{x — z) 

r+oo 

dzz^" / dx e^^^ dx u{x — z) 

J z 

r+oo 

dzz^°' / dxe^^^u{x~z) 



r(i-«) Jo 
r(i - a 

provided that (4.5) holds true. Furthermore, 
Til -a) Jo 

1*00 /' + C30 

dx e^^^uix) 



a«M(0=— ^/ dzz-"' I dxe'^^'uix- z) 



r(i-a) Jo 



iz z e 
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(4.6) ={-^irm^\Te^'"^^^m, aG(0,l) 

from the fact that formula (2.20) holds and provided that u E _L^(M-|_). From (4.G), 
the equation 

dug _ _ 
dt ~ ^''^ 

can be rewritten by passing to the Fourier transforms as 
which, in turns, leads to 

(4.7) ^{^, t) = i2J(0 exp (-i|^|"e-*^ ^) . 

If .4 = in (4.2), then from the operational rule (3.9), we obtain that 

u{x, t) — e^*^°=uo(a;) — uq{x — t) 

and thus 

poo 

(4.8) Ua(x,t) — e"*''^''^ uo{x) — I ds ha{s,t) uo{x — s). 

Jo 

The Fourier transform of the convolution (4.8) is therefore written as 

/•OO 

(4.9) ^(^,t)=/ ds/i„(s,i)e^«^^?(e) =^(OC(^,i) 



where hg{£^,t) = Eexpi^i3" is that in (2.12) and thus (4.9) coincides with (4.7). 
From this we obtain the claimed result. □ 

Remark 3. If we assume that u2{£,,Q) — 1, that is Ua{x,0) — S{x), than we get 
that Ua = ha which is the law of a stable subordinator Sjf, t > with a e (0, 1). 

Lemma 2. Let u : D C i—^-Rbea function such that u € Dom{a ■ V). For the 
fractional operator (3.11) we have that 

(4.10) (a'^^"M(^) = (-i a • I)" u(|). 



Proof. Let us consider the function u : D C M" M with 



u(x) 



0. 



We have that 



(a- V)"?/^) = / dxe'« 



: (a-V) 

r(i-a) 



/ dxe'«- rdzz-"e-^(-^)u(x). 
r(l - a) Jo 



From the Foubini's theorem we can write 



(eTvy^uiO = ^ f dzz-"" f dxe'«-^e-^('^-^)w(x) 

r(i - a) Jo Jo 
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where 

/ dx e*^'^e~^'-'*'^^u(x) = / dxi ■ ■ ■ dxn g^^^u{xi — aiz, . . . ,Xn — dnZ 

Jd JD+az 

^g«z(ai6+...+a„e„) f dxi---dXne'^-Mxi,---,Xn 

Jn 

From (2.20), by considering that 



we obtain 

r(l-a) Jo 
and formula (4.10) is proved. 

We present the main result of this section. 

Theorem 1. We have that: 

i) for /3 £ (0, 1] and a € (0, 1), the solution to the fractional equation 

(4.11) j^^ + (a.V)"^7/„,^(x,t)=0, xeM!;:,t>0 
subject to the initial and boundary conditions 

"a,^j(x, 0) = ■Uo(x), Ua,l3{O,t)^0 

is written as 

(4.12) u„,^(x,t) =Z^^(a-x,i), 

where is the solution to (2.2). 

ii) for /3 e (0, 1), the solution to the fractional equation 

(4.13) (af + (a-V))w,3(x,t) = 0, xeM';, i >0, ;3e (0,1) 
subject to the initial and boundary conditions 

z;/3(x,0) =(5(x), vp{0,t)^tl, u> -1 

is given by 

(4.14) v^{^,t)=f^W^^.^+^[-^). 

iii) for V = — /3, the problem ii) can be rewritten as 

(4.15) (^^ + (a.V)^«^(x,t)=0, xeM!^, i >0, /3e (0,1) 

subject to the initial condition 

«^(x,0) = ,5(x) 
and the solution (4.14) takes the form 

(4.16) w^(a • X, <) = ^^(a • X, t) 

where lp{x,t), {x,t) G (0, +oo)^ is the law of £^ , t>0. 
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Proof. Proof of i). From the previous result the Fourier transform of the equation 
(4.11) is written as 



The Mittag-LefHer function is an eigenfunction of the Dzerbayshan-Caputo frac- 
tional derivative and therefore we get that 

(4.17) iQ^it t) = Mi) Ep {-t^i-za ■ . 

This is equivalent to write Ua,i3 as follows 

u^,p{x,t)^ I ds?^(s,i)e-^("-^)°uo(x) 

dsU^{s,t)e-"'-^uoix), a,/3€(0,l) 



where, we recall that, 

/■OO 

(4.18) U^{x,t)= ha,{x,s)lfi{s,t)ds 

Jo 

is that in (1.4). From (3.9) we get that 

/■OO 



/•OO 

= / dsUf^{s,t)uo{xi - Sai, . . . ,an - QnSn) 

Jo 



and thus 

ill>(|,0= / dsW^(s,t)e'^=^-«^(a,...,en) 
Jo 

From (2.10), (2.11) and (4.18) we have that 

poo 

U^{-f,t) = I e'''''U^{x,t)dx 

J 

=i?0(-t^(-*7)"), 7GK 







and thus 



Proof of ii). The Laplace transform of vp is written as 

/•OO 

vp{-K,X)^ / e~^*u/3(x, t) 



dt 

dt 



k\T{-pk + v +1) 
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oo ^ 

=A— ^5:(-A^a.x)^-i 

fc=0 

(4.19) =A-''-iexp(-A^a-x) 
where, from (2.21), 

poo 

Jq 

We immediately get that 



vpix, A) = -A'' aj u^(x, A) 



and thus 

n 

dx 



II Q ft 



being ||a|l = 1. This means that 



(a • V)5^(x, A) = - A-^ ^^^(x, A) = ^ e^'^ (-9f Vf, (x, t)) 



where in the last identity we have considered the Laplace transform of the Riemann- 
Liouville derivative ([27]). 

Proof of iii). The Fourier transform of (4.15) is written as 

and leads to 

which coincides with (4.17) for a ~ 1. We recall that 

r ^"^^^"'^-^--^ (- J) = ^^-^^-^ (-^) ' ^ > 

where (see Beghin and Orsingher [5]) 



and, from (3.9), 



/ dxe'^ '^e '''^'^'^''vp{x,0) = dx e'^"t>/3(a;i - sai, . . . , a;„ - sa„, 0) 



where vpH) = 1- Thus, we obtain that 



dxe^«"'e-^(-^)z;^(x,0) J t-^W.p,,-^ (-|^) 



OO 

ds 
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= [from the operational rule (3.9)] 



where w^(x, 0) = (5(x) and x = sa if and only if x • a = s||a|p with ||a|p = 1. We 
obtain that 

{it^ ^■^)=£^ e'^ " ds r^W.p,i.p (-|^) (5(a • X - s)^ dx 

e*^'''w^(x,t) dx 



+ 



which concludes the proof. □ 
Remark 4. We remark that, for a — /3, formula (4.12) becomes 

sin^TT {a-x.f-'^tl^ 
up.,p[^,t) - ^^(j^.^)2/^^2(a.x)/3t/3cos/37r + t2/3' P ^ 

The law up ij{x,t)^ a; > 0, i > can be interpreted as the density law of the ratio 
of two independent stable subordinators both of order [3 G (Oil); see D'Ovidio 
[9]; James [16]. This ratio (formula (1.5)) does not depend on t > 0. 

Remark 5. For a ^ I and (3 G (0,1), the solution (4.14) (and (4.12)) with 

vpiQ^t) — t^^, can be written as 

vp{x, t) = t-Pw.p,i.p ("1^ ) =^P^'^' ^) 
where Ip is the law of £f , t > 0. 

Remark 6. For (3 — 1 and a £ (0, 1), the solution (4.12) takes the form 

(4.20) Ua (x, t) = Uq^i (x, t) = /iq (a • x, t) 
where ha is the density law of and solves 

(4.21) (^^ + (a-V)"^w„(x,i) = 0, (x,i) eM![ X (0,+oo) 
subject to the initial condition 

n 

(4.22) w„(x, 0) = 'uo(x) = Y[ M^k) 

fc=i 

with uo(xfe) — S{xk) for all A; 1, 2, . . . , n. Indeed, we have that 

/"OO 

(4.23) Ma(x,t) = e-*('"-^)°uo(x) = / ds ha{s,t) e''^''-^^uo{x:) 



and 

(4.24) ^;(^,i) =^(OC(a-|,i) 
where 

n 

(4.25) S5(0 = n^(^fc) 

fc=i 
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From (3.9) and the fact that uo{xk) — S{xk) for all k, we rewrite (4.23) as follows 



Ua(x, t) 



POO 

/ ds ha{s,t) S{x ~ sa) 
Jq 



By considering that x = sa if and only if a • x = s we obtain (4.20). From (3.9) 
and (4.22) we have that 

n 

fe=l 

and therefore we obtain the Fourier transform 

POO " 



a{i,t)^ e"^'^ I dsha{s,t)W_'^o{xk- sak)dyi 
fc=l -^0 

Formula (2.12) says that 

ha{ai-^,t)=cxp (-ila-ll^e-'^^THT) = exp {-t{-ia ■$,)") . 
By taking into consideration (4.10), from the fact that 

we obtain (4.21). 

From the previous result we arrive at the following statement. 
Corollary 1. For j3 G (0,1), the solution to the time- fractional equation 

(4.26) (af + (a- V))u/j(x,t) = 0, (x, t) e E'^ x (0, +oo) 
subject to the initial and boundary conditions 

i;^(x,0) = <5(x), V0{O,t)^tl^^ 

or equivalently 

(4.27) (^_ + (a.V)jz;^(x,t) = 0, (x, t) e M'; x (0, +^) 
subject to the initial conditions 

z;^(x,0) = 5(x), ^(x,0) = 0, < fc < Ml - 1 
(\-'\ is the smallest following integer) is given by 

(4.28) ^;^(x,t) = ^T4^_0,i_„^(-^). 

Proof. It suffices to consider v = — n/3 in Theorem 1 and the fact that 

|^(A) = A"^^I(A)- A-^-'^-^l^r)- 



□ 



FRACTIONAL ADVECTION DIFFUSION EQUATIONS 



19 



Remark 7. We point out that the normahzation of (4.28) given by 



(4.29) 



where a, 



p^(x,i;n) = ^W^, 



a • X 



X e M", i > 



in) 



ai ■ a2 ■ ■ ■ ttn, is a probabihty distribution whose one-dimensional 



marginals coincide with Ip, that is the density law of the inverse process £f , i > 0. 
Indeed, for the law pp{x, t; n), /3 e (0, 1), (x, t) e M" x M" , we can show that 



(4.30) 

and also that 
(4.31) 



/ P/3(x,t;n)dx = |5^(y,t;n- fc), ye 



t e 



/ p^(x,t;ri)dx = 1. 



To this end, we recall that (see for example Lebedev [19]) 



r(— a) 27ri Jq 



and we write 

p/3(x,t;n) 



l(n) 



(TI 



oo ^ / n \ 

.j=i ''j ) fc=i y j=i '^j I 



T{-^k + 1 - n/3) 



a 



(n) 



(n;uif)2-^c 

a(n) 1 
^(n"=iif)2^ 



k>0 



i=i j 

n 



3 = 1 ^3 



ds. 



By observing that 

we get that 

/ Pi3{x,t;n)dx 
Jo 



( 



4/9 



da;,- = — exp 



./3 



3-^3 

\ 3*3 / 



a(n) 1 1 
(n"=iif)^27*^ 



EUjJ 



1 j=i J I 

\ I 



I 



1 X ^ 1 X ^ ajXj 1 f 

k>0 1 = 1 ''3 ^ 



1 3=1 3 I 
\ 3*1 I 
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3^1 



\ ^ \ ^ ^j'^j 

\ 3^1 / 



;a-W^-/3,i-(«-i) 



37^1 



r(-/?fc + 1 - (n - 

IjXj 



. J=l 3 I 

\ 3+1 I 



which is the marginal law of 

p,3(x,t;n), X e M", t G M" 

and this proves (4.30). In order to obtain (4.31) it suffices to consider (4.30) and 
the fact that integrates to unity. Indeed, we recall that 



j(x, t; 1) = /^(x, t) = 



is the probability density on [0, +oo) of the process £f , t > 0. 

Remark 8. The distribution (4.29) can be regarded as the law of a process which 
is the inverse to the sum of stable subordinators. In particular, let us consider the 
stable random sheet 

(iiof^,...,„jofJ, x= (xi,...,a;„) eM'; 

with independent stable subordinators j^xj j J = 1: 2, . . . , n, of order e (0, 1). We 
are interested in studying the inverse multi parameter process 

5r^(<) = (ixf,...,„xf), t>o 

of the linear combination 



law p 



with aj > 0, yj, in the sense that 

(4.32) Pr{^pit) > x} = P{iXf > xi, . . .„xf > xn} = ^{^^(x) < 0- 
We first obtain that 

(4.33) Pr{Jf}3 (x) e ds}/ds = hp{s,a- x) 
Indeed, from (2.10), we get that 

(4.34) Eexp {-^Jfpi^)) = exp (-^^(a • x)) = Eexp (-^.$3^,,) 
which means that 

jr^(x) iof.^, aeM!^, xeM!^. 

From (4.32) and (4.33), we find the law of as follows 



P/3(x,t) 



dxi ■ ■ ■ dxn 

(-9)" 
dxi ■ ■ ■ dxn 



P{,y^li{^) < t} 
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= /* TT^^^i^W e ds} 
Jo dxi ■ ■ ■ dx„ 

From (2.10) and the fact that 

nOO / \ 1 /"OO 1 

J e-^'i^J hp{s,a-^)dsjdt=jj e-^"/i^(s,a-x)ds = -e-^"^-^)^" 
we arrive at the Laplace transform 

^ e-^'pf}{^,t)dt ^ e-^' (^J^ q^~^}q^ hp{s, a ■ x) dt 

dxi • • • dxn A 

n 

By recalling that a(„) = ai • • • a„ and inverting the above Laplace transform, we 
get that 

which is the claimed result. 

We now extend the results presented in Theorem L Let us consider the solution 
to the time-fractional problem 

(4.35) (af + (a- V))ii;^(x,t) = 0, x e M!;:, i > 0, /3 G (0, 1) 
subject to the initial and boundary conditions 

il^(x,0) = (5(x), il^(0,i) ^ir"*^"'' «eN. 
From Theorem 1 we know that 

(4.36) H^(x,t) = r-"'3-iM/-0,„_„^ ("ir) • 

Theorem 2. For n, m G N, we have that 

(4.37) il^(x, t) * il^(y, t) = ii;+™(x + y, t) 

where * stands for the Laplace convolution overt. Furthermore, formula (4.36) can 
he written as 

(4.38) ii^x,t) = {ui*---*un){x,t) 
where 

(4.39) ^i^(a;„^) = (-^) , x>0,t>0, i^l,...,n. 

Proof. We first assume that formula (4.38) holds true and we write 

0(x,i) = (lii * • • • * ■ti„)(x,t). 

We have that 

^ /-oo />oo 

mt)= ■■■ e-«-</)(x,t)da---rfe« 

JO Jo 
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OO /'OO 



e ^i=i (ui * • • • * Un){xi, . . . , a;„, dxi • • • dxr, 

'0 Jo 
= {u^ ...,£,„, t) 



where, from (1.6) and (2.4), 



1 



(4.40) u,{^^^t)^ I e-^^'=^u,{x,,t)dx,^ -E, 







and thus 

(4.41) 0(eA)=/ e-^*0(et)dt = n7T7r— 7T- 

Now we consider the function 

^(x,t)=r""^-il^_^,„_„0(^). 

From (1.6) we have that 

(4.42) ^(x,A)==/ e-^V(x,t)di = A"'^^"e^^''('^-^). 



As we can immediately check 

(4.43) ' - ■ ^ ^ 



/ e~^-^(^(x,A)dx^ (c 



..,Li(e. + A^a.) A^^(6 + A^a,) 
coincides with (4.41) and therefore cj) = ip. This concludes the proof. □ 
Remark 9. For n = 2 wc have that 

(^aixi)' (-aza^z)^' , .^p^.p , 



i=0 j=0 

OO OO 



where 



Z.Z.,,r(i_;3*-/3)j!r(i-;3j-;3)* ""^^ ^> 

r(2 - 2/3 - /3i - /?j) 



T{wi + W2) 

is the Beta function. From the fact that (see Lebedev [19, p. 115]) 

11 
r(a) ^ 2^ 

we can write 

00 00 



ii^ (X, t) ^t'-^^ y y t^^i^Ifl t^^^iifl 1. f ,2p+Pi+p,-2^s 

1=0 j=o ■' ■'^ 

i-2p±_ f^^ {~aixis^/t^y {-a2X2S^/t^y 2P-2 s , 

^ 1 n, A — n -J 



1=0 j=0 
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2i:i 



2TTt Jr f-' fc! 



,1-2/3 {~iaixi+ a2X2)/t'^)'' 
fro k\T{~l3k + 2-2f3) 



which is in accord with (4.36). 

Remark 10. We can obtain the previous result by considering that 



ill{^,t)^ / ui{S,i,s)u2{£.2,t- s)ds 
Jo 



where 

Thus, we arrive at 



u , (e, ,t) = Ep (-e, tP/a,)/a, , J = 1, 2. 



aia2 
where 



f l*oo poo 

E0{-^2S^/a2)Ep{-^i{t-s)'^/ai)ds^ / 6-^1-^^^(71,72)^71^72 



S(7i 72) = V V (-672/a2)' (-gi7i/«i)^" r\wu-sVPds 

/!r(/;3 + i) j!r(j/3 + i)* ^('/^ + 1,^/^ + 1) 

^ - - (-672/^2)' (-a7i/«i)^' 

ifefro ^- ^- r(2 + /3(/ + j)) 

^yy (-672/a2)' (-emM)^' ^ ( s-^e^ {t/ s^^'^^^' ds 
^ ^ l\ jl 2m Jr 

(-6/a272t^/s^)' (-a/«i7ii^/s_^ 



y" exp - ^(671/ai + 672/02)^ ds. 



1=0 j=Q ■' 

'2m 



By observing that 

t 

2TTi 



s exp - -^{^iji/ai + ^2^2/0,2) j ds 
e-,-2exp f_ (ei7i/ai+ 672/^2) \ 
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where C' is similar to C (see [19]), we obtain that 



27^^ 7c' 1 + 1 + ^^ 

1 /■ ^sf ,2/3-2 1/^1 V«2 



By inverting the Laplace transform we get the Laplace convolution of two Mittag- 
Leffler functions . 

5. Second order operators 

In the previous sections we have studied the fractional power of the first order 
operator (1.1). Here we present some result on the squared power of (1.1) and the 
fractional power i? € (0,1) of the negative Laplace operator —A which has been 
first investigated by Bochner [7]; Feller [12] and many other authors so far. 

Let us introduce the fractional power of (1.1) formally written as 

(n n „2 

with ^ ^ h j n and acting on the space of twice differentiable functions 

defined on M". We present some results for a = 2. As pointed out by Meerschaert 
et al. [21, 22], for a = 2, formula (3.8) leads to the directional derivative 

n n Ffi f( \ 

(5.2) Di,/(x) = VAV/(x) = EE«^.^-a^ 

i—l j — 1 * ^ 

where A = {aij}i<ij<n and 

e,ejM{de). 

=1 

Theorem 3. The solution to 

(5.3) (^^-(a.V)2^g(x,t)=0, (x, t) G M" x (0, +(X3) 
subject to the initial condition (7(x,0) — (5(x) is given by 

(5.4) g(x,i) = -l-W^_£._. ^ 



tl3/2 "-t,l-S \^ tl3/2 

Furthermore, formula (5.4) can be written as 



(5.5) gix,t) = I ^ /I l^{s,t)ds 



(2 



/47rs 



where 
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Proof. We recall that 

/.no 



From the fact that (see Lebedev [19]) 

1 f /.oo 



where K^, is the modified Bessel function of imaginary argument and ([19]) 

TT 



we obtain that 

(a>c) 

-At 



C30 fOO J— 

' e 1= 



e^-"' / lp{s,t)dsdt 

Jo Vins 



(|a.x|A^ 



0? V 4A'9 



4A/3 y V 2 

^ A 2 ^ exn I — la • xi; ~ 
which coincides with 



a-x|A2] ^ exp ( — [a • x|A' 



:^A2 exp ^— [a • x|A- 



oo 



l^(-|a.x|)'= p e-^H-P'-"-'^ dt 



Al-ig (~|a.x|A-^^- 

fc=0 



fc! 



-A 2 exp ( -|a • x|A = 



2 

This prove that formula (5.4) can be rewritten as (5.5). Formula (5.5) can be 
rewritten as follows 

/.oo / -1 /.+00 \ 

5(x,<)=y^ i—j e'^^<^'^--^e->'^dj\ lpis,t)ds 

lj3{s, t)ds ) dj 



1 

'2n 



-\-oo 

i7(a-x) / / ^-7^s; 



e 

— oo 

+00 



(5.6) =- / e-'^(--)£;4-tV)rf7 



OO 



/ ,-..(a.x)^^ (_,,^2) d7 = ^ / e--(--)(-7^)i?, (-t^^) d7 
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which coincides with 

^+00 



by considering that ||a|| = 1 and provided that "f'^Ep (— e V'ld'y). By taking 
the Fourier transform of (5.6) we obtain 

/• r+co 

git t) = / e'«-"g(x, t)db,= / 5{i ~ ^a.)Ep {~t^j^) dj 
where ^ = 7a if and only ifa-^ = 7||a|p=7 and therefore we get 

/+00 
<5(a . I - j)Ep i-tf'j') d7 = Ep (-t^(a • 1)^) . 
-00 

From the fact that 

/+00 r-\-oo 
e'^-''{sL-Vfu{x)dx = {-ia-£,f / e'^'^'uix) dx = - [a ■ £,f u{i) 
-00 J —00 

u G Dom (a ■ V), we get that the equation (5.3) becomes 

|^(^,t) = -(a.O' g{i,t) 

which leads to 

g{tt)^Ep(-tP (a.^^) 

by taking into account the initial condition (^(x, 0) = '5(x). We recall that the 
Mittag-LefHer function is an eigenfunction for the D-C fractional derivative ^ 
(see formula (2.5)). This concludes the proof. □ 

Remark 11. We notice that 

and therefore we can write the function (5.4) as follows 

g{x,t)= li{\a.- x.\,s)lp{s,t)ds. 
Jo 

Furthermore, the one-dimensional version of (5.3) is the fractional diffusion equa- 
tion 

g{x,t) = 0, xeR, t>0 



dtP dx^ 

subject to the initial condition 17(0;, 0) — 6{x) whose stochastic solution is repre- 
sented by the subordinated Brownian motion 

i?(£f), t>0 

as shown by Orsingher and Beghin [25]; Meerschaert et al. [24]. Furthermore, we 

notice that g(x,t) — h(\a ■ x\,t) and solves (5.3) whereas, l§_{a ■ x,t) solves 
2 2 

(5.7) [ -^-(a-V)Vii(a-x,t) = 0, (x, t) G K" x (0, +00) 

\dd J ^ 

as shown in Theorem 1. 
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Let us consider u e ^ where .5^ is the Schwartz space of rapidly decaying C°° 
functions in M". The fractional power d G (0, 1) of the negative Laplace operator 
— A is defined as follows 

C f u(x + y) + m(x - y) - 2w(x) 



2 7k. |y|"+2'' 

where C is a constant depending on (i?, n) and, "p.w." stands for "principal value". 
An alternative way (in the space of Fourier transforms) to define the fractional 
power of —A is 

-(-A)''u(x) = ^ / e-^^-Uruii)dC 
Indeed, we have that 

(5.8) / e'«-(-A)''u(x)(ix=||^f^2(0. 
The solution to the fractional equation 

(5.9) |^(x,i) = -(-A)''w(x,i), xeM", t>0 
in the Sobolev space 

(5.10) i?om(-(-A)^) = (uei^^M") : / (1 + |l||l2*)|u(|)|2d| < oo 



and subject to the initial condition uo(x) = '5(x), represents the law of an isotropic 
]R"-valued stable process, say S2i?(i), i > 0, with characteristic function 

(5.11) Ee'«-^^^W =e-*ll«ll''. 

Thus, from our viewpoint, the isotropic stable process (i), i > 0, G (0, 1) is 
the stochastic solution to (5.9) subject to the initial datum mq = b. Obviously, for 
1? = 1, we obtain S2(i) = B(i), i > which is the n-dimensional Brownian motion. 
We now present the main result of this section. 

Theorem 4. For -d e (0, 1), a,/3 S (0, 1), a e M" such that \\a.\\ = I, the solution 
w = w(x, t) to the fractional equation 

(5.12) ^^_(_A)''u;-(a.V)"^ 

for (x,i) G M" X (0, +cxd), subject to the initial condition wq — S, is given by 

/ dz dsT2i){x ~ s, z) ha{a ■ s, z) lfj{z,t) 
Jo Jw'i 

where 72i?(x,i), x £ M", t > is the law of the isotropic stable Levy process S2^{t), 
t > 0; ha{x,t), X e [0, +00), t> Q, is the law of the stable subordinator Sjf , t > 0; 
ll3{x,t), X £ (0, +00), t > 0, is the law of the inverse process £f , t > 0. 

Proof. First, from formulae (3.11) and (3.9), we evaluate the Fourier transform of 
a function u G L^(M") given by 

f e*«-^(a- V)"it(x)dx= / e»$-x ' ^) f ds s-"e-'^''-'^\{x)dx 
7r" 7r" r(l - a) Jo 



(5.13) w(x,t) 
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Ha-O 



r(i- 


a) 


(— ia 





r(i- 


a) 


(~ia 





r(i- 


a) 



dss'" 



(5.14) 



= (-za.e)"2(0 



(where we used (2.21)) which is in accord with (4.10) for functions in L^{MJ^). From 
(5.8) and (5.14), by passing to the Fourier transform, the equation (5.12) takes the 
form 



(5.15) 



w , 



From the fact that 



/ e'« ''r2i^(x-s,z)dx = e'«-" / e'« ''r2^(x, z) dx e*«-" e-"ll«ll 



we get that 



w(^,t) = j e^«"u;(x,t)dx 



From Remark 6 (formula (4.24) in particular) we have that 

Jq 

and, from (2.11), 

(5.16) Witt) = Ep {-t^Wif"^ ~ t^i-ts. ■ 0") . 

From the fact that 



dtP 



{-t\) = -CEp{-t^O, C>o 



/5/ 



which means that the Mittag-Lefher is an eigenfunction for the Dzerbayshan- 
Caputo time-fractional derivative, we arrive at 



dtp 



(tt) = -(llllP'' + (-*a-0") wii,t) 



which coincides with (5.15). This concludes the proof. 



□ 



Remark 12. The distribution (5.13) can be regarded as the law of a subordinated 
E"-valued stable process with stable subordinated drift given by 

W(t) = S2^(£f)+ajo2,, t>Q. 



The characteristic function is given by 



EgiC-S2<)(T)+j«.ai5? 
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=E 



e-Til€f''Ee-(-^€-a)i5? 



T = £f 



= (by (2.10)) 



=(by (2.11)) 

which coincides with (5.16). 

Remark 13. We remark that, for a = (3, the stochastic solution to (5.12) becomes 

S2^(£f) + ai4j' *>0. 

where the ratio of two independent stable subordinators j^^, j = 1,2, is indepen- 
dent from t and possesses distribution 



(5.17) Pr{iS)^/2Sj'^ edx}/dx = 



sin /Jtt 



x>0, t>0. 



TT x^l^ + 2xf COS /37r + 1 ■ 
This has been pointed out also in Remark 4. 

Remark 14. For ^? = 1 in (5.12), we immediately have that the fractional equation 



dtp 



has a solution which is the law of the process 

B(£f)+ai5^,, t>0 
(B is a Brownian motion) and can be written as follows 



ha{SL- s,z)lp{z,t). 



w{x,t) ~ dz ds — 

Jo Jr^ V 47rz 

For a — >■ 1 we have that Sjf — >■ t which is the elementary subordinator (obviously 
ha{x,t) — >■ S{x — s)) and therefore the solution to 



dt0 



Aw — (a • V)w 



is written as 



w{x, t) 



dz I ds 



ds 



where we used the fact that a • s = z G M+ iff s = az or, equivalently 



w{ii, t) ^ dz 
Jo 



ds 





e 






J Atkz 












/47rz 




_ iff s 










^/Attz 



(5(a • s — z) /^(z, £) 
8(s — az) Z^(z, €) 



(5(a • s — z) Z^(z, i) 
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r g 4(as) 

■ / ds dp{a-s,t). 
Jb.1 V47r(a-s 



v/47r(a • s) 

Finally, for d = 1 and a = 1, we get that 



°° e- 



w(x,t) = / ^ — li3iz,t) 

V47rz 





is the law of the n-dimensional subordinated Brownian motion with subordinated 
drift 

B(£f)+a£f, t>0. 
If /3 — >■ 1, then £f — >• i and lp{x, t) — >• 5{x, t). The solution to 

dw 



„ — Alt; — (a • V)w 



is given by 



which is the law of 



zi;(x, t) 



/A-Kt 



B(t)+at, t>0. 



Remark 15. For the sake of completeness we also observe that, for a = 1, equation 
(5.12) becomes 

(5.18) _^_(^A)''^-(a.V)z. 

whose stochastic solution is given by the subordinated stable process with drift 

(5.19) S2i,(£f) +a£f, t > 0. 
Indeed, for a — > 1, we have that j^" — > t. 

Remark 16. We remark that ([28; 29]) the transport equation 

du 

= Au — Am + XKu 

at 

where 

^ 

k=l " 

and K is the Frobenius-Perron operator corresponding to the transformation T(x) = 
X + f{x) has a solution which is the law of the solution to the Poisson (Nj) driven 
stochastic differential equation 

(5.20) dXt = a{Xt)dt + /(Xt)dNt. 

Remark 17. We recall that (see for example Beghin and Orsingher [5]) the frac- 
tional Poisson process N'f = 7V(£f ) is a renewal process with 

Pr{A/-f = fc}=pf(i)=Ep,(£f) 

where Pk{t) is the distribution of the Poisson process N{t), t > and probability 
generating function written as 

(5.21) Ez< = Efi (-A(l - z)t^) . 
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From (5.21) and (2.5) we can write 

pip , » 

Furthermore, we recall that 

(5.23) Ez^(*) = e-^*(i-^) 
and 

(5.24) Ee-*«^(*) = exp (-^(1 - £-'«)) . 
Theorem 5. The stochastic solution to 

(5.25) ^(x,i) = -((a.V)" + A(/-e-i-^) jw(x,t), (x, t) e M'; x (0, +oo) 

where a, /3 E (0, 1] and e^"'"'^w(x, t) = w(x— 1, t) is the shift operator, is the process 

(5.26) Yt =N(£f)+aio^,, t > 0. 

Proof. We have that 

Ee*«^* =Eexp (-A(l - e*«)£f + ia • 

=Eexp (-A(l - e'«)£f - (-za • 0" ) 
=£;4~A(l-e*€)t^~(~*a.0"i^) 

is the characteristic function of (5.26). From (2.5) we obtain that 

(e t) = (-A(l - e^«) - (-^a • 0") w(|, t) 
where, from Lemma 2, (a • V)"w(^,t) = (— ia • ^)"w(^,t) and 

/" e''«''A(/-e-i-^)w(x,t)rfx=A(l- /" e^«''w(x - 1, t) dx 

=A(l-e*«)w(^,t). 

This shows that (5.26) is the solution to the Poisson driven stochastic differential 
equation whose density law solves (5.25). □ 

Remark 18. Let us consider the equation 

(5.27) -g^ix,t)^~[d^ + X^ ^)w(x,t), (a;,i) eM+ X (0,+oo) 

which coincides, for t — I, with (5.25) in the one-dimensional case and can be 
rewritten, for r ^ 0, as follows 

, A 



(5.28) [W^^") = "7 " " ^' • 
Formula (5.28) is the governing equation of the one-dimensional process 

(5.29) rt = TiV(T-i£f)+io^,, t>0. 
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For r — >• 0, we obtain that 

(5.30) (^^+9^^w(x,t) = -A9,w(x,t) 
is the governing equation of 

(5.31) Yt = \t + S^%,, i>0 
which becomes, for a = /3 G (0, 1), 

(5.32) Yt=(^\+ i > 

where the ratio of independent stable subordinators has distribution (5.17). 
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